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This is a great website: Images of Mathematicians on postage stamps 

http://ieff560.tripod.com/stamps.html 

with not only the images but much information about mathematicians from all over the world 


Download two Computer Algebra Systems free of charge: 

Octave and Xoctave: http://www.qnu.orq/software/octave/ 

XCAS: 

Of Aide for free software XCAS , aid written by Renee de GRAEVE Universite J Fourier 
GRENOBLE FRANCE. This AIDE is enormous. I have not the courage for translate it for you. 
You may find it on the site 

www-fourier.uif-qrenoble.fr/~parisse/qiac fr.html 

Thanks to Robert Setif for this information 


A wonderful book about Strange Attractors can be downloaded free of charge (together 
with programs) from 

http://sprott.physics.wisc.edu/sa.htm 

http://sprott.physics.wisc.edu/fractals.htm 

I plan an extended DNL-paper on this fascinating topic. “Create your own Strange Attractor /” 
(See two graphs on page 43.) 

This is a website on Curves (in German) 

http://www.blaesius-anne.de/arbeit kurven/kurvenGla.html 

SOME RELATED PUBLICATIONS FROM PREVIOUS FORMATEX ACTIVITIES -Book: 

“ Research , Reflections and Innovations in Integrating ICT in Education 

http://www.formatex.orq/micte2009/volume1.htm 

http://www.formatex.orq/micte2009/volume2.htm 

http://www.formatex.orq/micte2009/volume3.htm 

Zwei sehr schone Bucher: 

Rainer Muller, Klassische Mechanik, De Gruyter 2010, ISBN 978 3 1 1 025002 2 

Georg Glaeser, Wie aus der Zahl ein Zebra wurde, ein mathematisches 
Fotoshooting, Spektrum 2010, ISBN 978 3 8274 2502 7 

For TI-Nspire novices: 

Steve Oullette, TI-Nspire for Dummies, Wiley Publishing 2009, ISBN 978 0 470 37934 9 

Written by a math teacher, TI-Nspire For Dummies helps students and educators alike take full 
advantage of this outstanding tool for teaching and learning math. 
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Liebe bUG-Mitglieder, 
und wieder war ich nicht rechtzeitig mit 
dem bNL! Wir haben in den 50 Seiten 
zwar „nur" drei langere Artikel, aber diese 
haben in einem regen und langen Aus- 
tausch von emails eine bemerkenswerte 
Eigendyamik entwickelt. Viele Fragen, 
Ideen, Anderungen und Verbesserungs- 
vorschlage gingen hin und her. Wir hoffen, 
dass sich die Bemuhungen fur die Leser 
der Beitrtige lohnen. 

buncan Me Dougall zeigt interessante 
Zusammenhange in Polynomen mit ganzen 
Koeffizienten und rationalen Extremwer- 
ten. Piotr Trebisz nutzt die ganze „power“ 
eines CAS fur die Berechnung und bar- 
stellung von Schneckenhausern. Robert 
Setif stellte mir ein nettes Problem, das 
ich dann doch in seinem Sinn gelost haben 
diirfte. 

bietmar Oertel hat zwei inhaltsschwere 
Papiere geschickt, die demnachst im bNL 
zu finden sein werden. Roger Folsom stellt 
im User Forum eine Frage an die bERIVE 
User, ein longer Brief von ihm steht fur 
das ntichste User Forum bereit. Auch mit 
bavid Sjostrand hat es einen ausgedehn- 
ten - mdglicherweise noch nicht beende- 
ten - Briefwechsel bezuglich des „Bewei- 
ses" aus dem bNL#80 gegeben (Seite 
40). 

Ich mbchte Euch besonders auf die In- 
ternet Links hinweisen. Zwei mtichtige 
CAS sind frei verfugbar. bas Buch von 
Julien C. Sprott uber „Strange Attrac- 
tors" ist eine Schatzkiste fur "Frakta- 
listen". bie deutschen Bucher auf der 
Info-seite sind ein wahrer Genuss fur 
Auge und Hirn. 

Leider gibt es dieses Mai nur wenig zum 
TT-Nspire. Version 3 kann beeits von der 
TT-homepage bezogen werden. bie neuen 
Moglichkeiten wollen wir im nachsten bNL 
besprechen. 

Ich wiinsche Euch eine schone Zeit und 
verbleibe mit den besten GriiBen 


bear bUG Members, 

and again I was not in time with the bNL! 
We have 50 pages and „only" three more 
extended articles. But these developped a 
remarkable momentum from an intense 
exchange of emails. Many questions, ideas, 
changes and proposals for improvements 
went to and fro. We hope that our ef- 
forts will be appreciated by the readers. 

buncan Me bougall shows interesting rela- 
tionships in polynomials with integer coef- 
ficients and rational turning points. Piotr 
Trebisz uses the full „power" of a CAS 
for the calculation and presentation of 
snail shells. Robert Setif presents a nice 
problem, which I could solve in his sense. 

bietmar Oertel sent two very contentful 
papers which can be found in the next 
newsletters. Roger Folsom presents a 
request in the User Forum, a long letter is 
ready for publication in the next User 
Forum. I had an extended - possibly not 
yet finished - email exchange with bavid 
Sjostrand concerning the „Proof" from 
bNL#80 (page 40). 

I'd like to point especially to the internet 
links. Two powerful CAS can be 
downloaded free of charge. Julien C. 
Sprott's book about "Strange Attractors" 
is a treasure box for "fractal ists". The 
German books presented on the info page 
are a true pleasure for eyes and brain. 

Unfortunately we have little stuff for the 
TT-Nspire this time. Version 3 can be 
downloaded from the Tl-website. We will 
talk about the new features in the next 
bNL. 

I wish you a good time and remain with my 
best regards, 



Download all ZWL-Derive- and Tl-flles from 

http : / / www . austromath . at/ dug/ 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a contents of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
77-CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bdhm 
D'Lust 1, A-3042 Wurmla 
Austria 

Phone: ++43-(0)660 3136365 

e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: June 2011 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Tools for 3D-Problems, P. Liike-Rosendahl, GER 

Financial Mathematics 4, M. R. Phillips 

Hill-Encription, J. Bohm 

Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 
An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 
Overcoming Branch & Bound by Simulation, J. Bohm, AUT 
Graphics World, Currency Change, P. Charland, CAN 
Cubics, Quartics - Interesting features, T. Koller & J. Bohm 
Logos of Companies as an Inspiration for Math Teaching 
Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 
BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 
Truth Tables on the Tl, M. R. Phillips 

Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 
Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZ & Rob Gough, UK 

A Conics-Explorer, J. Bohm, AUT 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

Treating Differential Equations (M. Beaudin, G. Piccard, Ch. Trottier) 

Structured Combinatorics, D. Oertel, GER 
A new approach to Taylor Series, D. Oertel, GER 
Statistics with Tl-Nspire, G. Herweyers, BEL 
Cesar Multiplication, G. Schodl, AUT 
Find your very own Strange Attractor, J. Bohm, AUT 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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I invite you browsing in higher mathematics journals, “emis” opens a spectacular view to the work of 
many mathematicians free of charge. You can find a rich collection of websites hosting mathematics 
journals from universities and mathematical societies from all over the world. 


The Electronic Library of Mathematics 



http://www.emis.de/elibm/iournals/index.htnnl 


Acta Mathematica Academiae Paedagogicae Nyireghaziensis 
Acta Mathematica Universitatis Comeniana, 

Acta Universitatis Apulensis, Alba Julia, Romania 

An International Scientific Journal of Sapientia University of Transylania, Romania 
Analele Stiintifice ale Universitatii Ovidius Constanta, Romania 
Annales Academice Scientiarum Fennicce, Helsinki 
Annales Mathematicae et Informaticae, Esterhazy Karoly College Eger 
Applied Mathematics E-Notes, Taiwan 

Archivum Mathematicum, Faculty of Science ofMasaryk University Brno 

Advances in Geometry 

Algebra Montpellier Announcements 

Algebraic and Geometric Topology (see also: http://msp.warwick.ac.uk/gtp/) 

Balkan Journal of Geometry and Its Applications (Balkan Society of Geometers) 

Beitrage zur Algebra und Geometrie / Contributions to Algebra and Geometry 

Bulletin of the Venezuelan Mathematical Association 

Sociedad de Estadistica e Investigacion Operativa 

Bulletin of the Belgian Mathematical Society Simon Stevin 

Bulletin of the Malaysian Mathematical Sciences Society 
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Bulletin of the Malaysian Mathematical Sciences Society 

Bulletin, Classe des Sciences Mathematiques et Naturelles, Sciences mathematiques 
Divulgaciones Matematicas, University of Zulia, Venezuela 
Documenta Mathematica - Journal der deutschen Mathematikervereinigung 
Electronic Journal of Probability, University of Washington 

The Electronic Journal of the Argentine Society for Informatics and Operations Research 
The Electronic Journal of Combinatorics, University of Pennsylvania 
The Electronic Journal of Differential Equations, Southwest Texas State University 
The Electronic Journal of Linear Algebra, Israel 

Electronic Journal of Qualitative Theory of Differential Equations, Univ. of Szeged, Hungary 

Electronic Transactions on Numerical Analysis, University of Kent, OH, USA 

Ensaios Matematicas, Brazilian Mathematical Society 

Hindawi Open Access Journals, A bundle of various journals 

International Journal of Open Problems in Computer Science and Mathematics, 

INTEGERS - Electronic Journal of Combinatorial Number Theory, Univ. ofWest-Georgia 
Journal of Applied Analysis, Lemgo, Germany 
Journal of Complex Analysis, Lemgo, Germany 
Electronic Journal for History of Probability and Statistics, 

Journal for Geometry and Graphics, Lemgo, Germany 
Journal of Complex Analysis, Brown University, USA 

Journal of Inequalities in Pure and Applied Mathematics, Victoria University of Technology 

Journal of Integer Sequences, University of Waterloo, Canada 

Journal de Theorie des Nombres de Bordeaux 

Living Reviews in Relativity, Max Planck Institute 

Lobachevskii Journal of Mathematics, Kazan State University 

Matematicki Vesnik, Mathematical Society of Serbia 


Mathematica Bohemica 


go to page 50 
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Diophantine Polynomials (Part 1) 

Duncan E. McDougall, duncanemcdougall@hotmail.com 

a) 

What is a Diophantine Polynomial? It is a polynomial of degree 2, 3 or 4 which is factorable in 
the set of integers and whose derivative is factorable in the set of rational numbers. Why we want to 
discuss them is to facilitate curve sketching. 

The polynomials which we are about to examine can be used for both the grade 1 1 and the cal- 
culus student, because the intercepts are easy to find, and the y-values for the maxima and minima are 
shared among the families of curves. For example, we can ask a grade 11 student to sketch 
y = x 3 + x 2 - 16x - 16 by finding both the x and y intercepts. We can use the very same polynomial 
for the calculus 12 student who can not only find the intercepts easily, but can more readily find the 
x- values for both maxima and minima because the derivative is easy to factor. 

My belief is that students should learn a complicated algorithm in simple progressive steps us- 
ing straightforward numbers. Diophantus worked with integers and rational numbers only. Pedagogi- 
cally, Diophantus was really on to something because he created methods which involved a lot of 
processing and sequencing while focusing on whole numbers. The distractions I refer to in curve 
sketching are complex and irrational numbers. It is difficult enough to learn some five to eight steps 
gathering enough data in order to effectively and accurately sketch a cubic, quartic, or quintic poly- 
nomial and/or a rational expression which may involve a diagonal asymptote without difficult-to- 
work-with numbers. If the student has the burden (when first learning the process) of working with 
irrational or complex numbers, along with concentrating on behavior of the curve and concavity, then 
he/she might simply declare 


( 2 ) 

“whatever” and drop the task. If the numbers are whole or integral (Diophantus), then his/her 
focus remains where it should be- on the algorithm. The task of the educator is to demonstrate algo- 
rithms in such a way that the student can master the process in sequence. The solution is to stick to 
the Diophantine process and to model examples that involve process and sequencing without getting 
tangled up with irrational numbers. To some readers this may be a no-brainer, but it is not as simple 
as it sounds to find cubics, or quartics with single integral roots whose derivatives have single rational 
roots. Finding these involved testing hundreds of polynomials using DERIVE, as I was determined to 
find easy-to-calculate polynomials, which would facilitate graphing curves like y = x 3 + 1 lx 2 + 24x 
without worrying about irrational and complex numbers. There was another challenge of course, and 
that was to keep the constant of the polynomial relatively small so that working without a calculator 
would not be arduous. 

Another aspect of this approach with whole numbers is that when the student knows that the 
numbers are designed to work, learning of the method or algorithm remains the priority. The student 
also knows that there is something wrong if the numbers do not work. It’s kind of a security blanket 
for the beginner, but it eliminates doubt that so often takes away confidence in ability and perform- 
ance. Later on, after mastering the technique, the student gains confidence through the ease of this, 
and therefore he can tackle problems with both irrational and complex numbers. 
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It is my objective to propose families of cubics and quartics which are factorable in the integers 
and whose derivatives are factorable in the set of rational numbers. I will also 


(3) 

propose methods using DERIVE by which you can construct your own polynomials. We’ll start with 
the very basic table of linear and quadratic polynomials then lead up to the cubics and quartics. I will 
end the paper with a brief discussion of the quintic, which should have worked but didn’t. 

The following diagram (Table 1) contains all the various linear and quadratic forms along with 
the general set of cubics. 


Using Table I 


For example, take a cubic of the form (x + a 2 )(x + 2ab}(x + h 1 j whose derivative has ra- 
tional roots. Choosing any integers a = 1 and b = 3 for example, our new polynomial is 
(x + l)(x + 6)(x + 9) with roots -1,-6, and-9. The differential form is 3x 2 + 32x + 69 whose 

23 
3 ' 


roots are -3 and 


(4) 


Table I 


Family Function 

Roots 

Derivative 

Roots 

Conditions on 
coefficients and 
constants to 

have integral 

roots 

a 

none 

0 

none 

not applicable 

a x 

x = 0 

a 

none 

not applicable 

ax + b 

_-b 

a 

a 

none 

not applicable 

( x + a) 2 

x--a 

2 (x + a) 

x--a 

none 

x + ax 

i 

o' 

ll 

X 

2 x + a 

a 

x = — 
2 

a must be even 

x 2 + x(a + b) + ab = 
= (x + a)(x + b) 

x = —a,—b 

2x + a + b 

-a-b 

x- 

2 

a and b are both 
odd or both even 

acx 2 + x(ad + be) + bd = 
= (ax + b)(cx + d) 

b d 
a c 

2 acx + ad + be 

-ad - be 

x = 

lac 

a #0, c 7^0 
ad + be must 
either equal ac 
or be an even 
multiple of it 

(x + a ) 3 

x--a 

3(x + a) 2 

x--a 

none 
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Family Function 


(. x + a) 2 (x + b ) 


x(x + a)(x + b) 


(x + a )(x + 2 ab)(x + b ) 


(x + l)(x - # )(x + a) 


Roots 

Derivative 

Roots 

x = -a, 

-b 

2(x + a) (3x +2 b + a) 

-2b -a 

x = - a , 

3 

x - 0, 

— a,—b 

3x 2 +2 x(a + b) + ab 

-( a + b)± ^ [a ‘ 
x - 

3 

x = - a 2 , 
-2 ab,-b 2 

3x 2 + 

+x(2a 2 + 4 ab + 2 b 2 ) + 
+ab(2a 2 +ab + 2b 2 ) 

x = - ab , 

-2a 2 - ab- 2 b 2 
3 

x = -1, 
a, -a 

3x 2 -2 x-a 2 

l + Vl + 3a 2 

x = 

3 


Conditions 
on coeffs & 
consts to 
have integral 
roots 

2 b + a is a 
multiple of 3 
a 2 - ab + b 2 
equals zero 
or a perfect 
square 

2a 2 -ab + 2b 2 
must be a 
multiple of 3 

4+3 a must 
be a perfect 
square (< a 
=0,1,4,15,...) 


The polynomials in Table II consist of the particular numerical families with single roots. 
These are the ones that are ready to use in your classroom today. 

As we observe the families in Table II (page 9), it is hard not to notice the pattern 8, 15, 21, 30, 
35, and 36. It is a quadratic arithmetic sequence, whose elements (except for a couple) all work as 
families of curves. 

In terms of methods for single roots let us begin by entering the form x(x + a)(x + Z?) into 

DERIVE. This guarantees a factorable form. Press C for Calculus and differentiate. The resulting 
form is put in function form as DECLARE. Now we can either guess values and hope that our quad- 
ratic is factorable, or fix a value for “a”, and then guess values for “b ” until the quadratic is factor- 
able. The question is, do we have anything to guide our guessing? In fact, we do. Visually, the val- 
ues of “x” for maxima and minima will occur between the first and last x-intercepts. Hence, if we 
were to choose 0 and 8 as two of our first and last roots, we would know that the third one must come 
between them. It’s just a question of leaving enough room between the roots so that the critical points 
can occur as integers and/or rational numbers. Algebraically, we enter x(x - <z)(x - 8) into DERIVE, 

and then differentiate giving 3x 2 + 2 x(a - 8) + 8 a. Since 


we have a quadratic; the discriminant B 2 - 4 AC must equal a perfect square in order to be factorable. 
Using the command DECLARE, we set / (a) = 4 a 1 - 32 a + 256 = 4 (a 2 - 8a + 64 ) and evaluate (or 

use the TI83 where second function gives “TABLE” and we search it for perfect squares). Both 3 and 
5 come up quickly implying both x(x — 3)(x — 8) and x(x — 5)(x — 8) have derivatives whose roots 


are rational. 
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#1: x- (x - a) ■ (x - S) 

d 

#2: — (x- (x - a) ■ (x - 8)) 

dx 

2 

#3: 3-x - 2-x-(a + S) + S-a 

2 

#4: FACTOR (3 -x - 2-x- (a + 8) + 8-a, Rational, x) 

2 

#B : 3-x - 2-x- (a + S) + S-a 

You may like to continue working in DERIVE; then use the TABLE-command in connection 
with SELECT in order to find (more) values for a leading to a perfect square, Josef: 

2 

#6: 50LVE(3-x - 2-x-(a + 8) + S-a, x) 

2 2 
-/(a - 8-a + 64) - a - S /(a - S-a + 64) + a + S 

#7 : x = - v x = 

3 3 

f 2 T 

#8: SELECT| FLOOR C/(v )) = ), v, TABLE(a - S-a + 64, a, -100, 100) 

J 

d 2 

#10: — (x-(x + 7) ■ (x - &)) = 3-x - 2-x - 56 

dx 

2 r 14 ^ 

#11: SOLVE (3-x - 2-x - 56, x) = x = v x = -4 

L 3 J 

L 15 169 J 

I don’t pretend to have all the families necessarily, but applying translations to any given family will 
yield many polynomials. The following is a small sample arrived at by adding a constant to all the 
terms: 

given family x(x-3)(x-8) 
add 1 (x + l)(x-2)(x-7) 

add 2 (x + 2)(x- l)(x- 6) 

add 3 (x + 3)x(x-5) 

add 4 (x + 4)(x + l)(x-4) 

add 5 (x + 5)(x + 2)(x — 3) 

add 6 (x + 6)(x + 3)(x — 2) 

add 7 (x + 7)(x + 4)(x-l) 

(x + 8)(x + 5)(x), etc. 


-7 169 

0 64 
3 49 

5 49 

S 64 


add 8 
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Interestingly enough, the entire above (all members of the famly) shares a maximum height of 
400 5 4 20 , 


27 3 3 3 


and minimum low of -36, and the difference between their corresponding 


x-coordinates is exactly — . (The students can prove this, Josef) 

A linear relationship exists between these values and those found in the quartics. We shall ex- 
plore this after exploring quartic family of curves. 


Family Function 

Roots 

Derivative 

Roots 

Transformation 

x(x + 3)(x + 8) 
x 3 + 1 lx 2 + 24x 

0, -3, -8 

(3x + 4)(x + 6) 
3x 2 + 22x + 24 

4~ 6 

(x ± k)(x ±3 a± k){x ±8 a + k) 

x(x + 5)(x + 8) 
x 3 +13x + 40x 

0, -5, -8 

(3x + 20)(x + 2) 
3x + 26x + 40 


(x ± k){x ±5 a± k){x ±8 a±k) 

x(x + 7)(x + 15) 
x 3 + 22x +105x 

o, -7, 

-15 

(3x + 35)(x + 3) 
3x 2 +44x + 105 


(x ± k)(x ±la± k){x ±15 a±k) 

x(x + 8)(x + 15) 
x 3 + 23x 2 + 120x 

0, -8, 

-15 

(3x + 10)(x + 12) 
3x 2 +46x + 120 


(x ± £)(x ± 8a ± k)(x ±15 a±k) 

x(x + 5)(x + 21) 
x 3 + 26x +105x 

0, -5, 
-21 

(3x + 7)(x + 1 5) 
3x +52x + 105 


(x ± k)(x ±5 a± k){x ±21 a + k) 

x(x + 16)(x + 21) 
x 3 +37x 2 +336x 

0,-16, 

-21 

(3x + 56)(x + 6) 
3x 2 +74x + 336 

56 * 
"T’" 6 

(x ± k){x ± 16a ± k){x ±21 a±k) 

x(x + 26 - a)(x + 26) 
a 3 + x 2 (52 - a) + 26x(26 - a) 

0, a - 26, 
-26 

3x 2 +52x + 2 6a -a 2 

not rational 


x(x + 14)(x + 30) 
x 3 + 44x 2 + 420x 

0,-14, 

-30 

(3x + 70)(x + 6) 
3x 2 + 88x + 420 

70 ^ , 
“T’“ 6 1 

(x ± k){x ± 14a ± k){x ± 30a ± k ) 

x(x + 16)(x + 30) 
x 3 + 46x 2 + 480x 

0,-16, 

-30 

(3x + 20)(x + 24) 
3x 2 +92x + 480 

20 „„ 
-y - 24 

(x ± k){x ± 16a ± k){x ±30 a±k) 

x(x + 33-a)(x + 33) 
x 3 +x 2 (66-a) + 33x(33-a) 

0, a - 33, 
-33 

3x^ + 2x(66 - a) + 
+1089 - 33a 

not rational 


x(x + ll)(x + 35) 
x 3 + 46x +385x 

0,-11, 

-35 

(3x + 77)(x + 5) 
3x +92x + 385 

-11-5 

3 ’ 

(x ± k){x ± 1 la ± k){x ± 35a ± k ) 

x(x + 24)(x + 35) 
x 3 + 59x 2 + 840x 

0, -24, 
-35 

(3x + 28)(x + 30) 
3x 2 + 1 1 8x + 840 

28 

,-30 

3 

(x ± k){x ± 24a ± k)(x ± 35a ± k) 

x(x + 36 - a)(x + 36) 
x 3 + x 2 (72 -a) + 36x(36 - a) 

0, a - 36, 
-36 

3x + 2x(72 — a) + 
+1296-3 6a 

not rational 
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Having fully explored the cubic, the quartic family of curves presented quite a challenge be- 
cause there would be three roots, (other than zero), to find. Visually, I opted for a span of 7, (one less 
than the 8 for cubics) entered x(x-a){x-b){x-l} into DERIVE, fixed a = 3 (only because it had 
worked with the cubic), took the derivative and evaluated “b ” from 1 to 7 hoping some value “b ” 
would work. The derived form was 4x 3 + x 2 (-3 b - 30) + x(20 b + 42) - 21 b and by declaring “f’ as the 
function I simply tested values for “b” and systematically factored (pressing F). To my great delight 
b — 4 worked, giving 2(x-6)(x-l)(2x-7). Observing 3 and 4 together, I acted on a hunch that 

Pythagorean Triples might work. So following in the footsteps of Diophantus, I tried triplets begin- 
ning with odd numbers and even numbers, and they worked beautifully. An added bonus were those 
triplets with consecutive legs such as 20 - 21 - 29 and 119 - 120 - 169, etc., which also worked won- 
derfully. The patterns appear in Table III. 

The numerical families of the form x (x + a) (x + b) (x + c) appear in Table IV. The numerical 
families for the form x 2 (x + a)(x + Z?) appear in Table V. 

Editor’s comment: How to use Table II for obtaining a bundle of cubics giving extremal values 
with rational x-coordinates: 

#1: p(x, a, k) := (x + k)-(x + 3- a + k)-(x + 8- a + k) 

Two examples using Table II 

' (x t 3)-(x - 3).(x + 13) ' 

(x - 2) ■ (x + 4) ■ (x + 14) 

(x - 1) ■ (x + S) ■ (x + IB) 

#2: VECT0R([p(x, 2 , k) ] , k , -3 , 3) = x- (x + 6) ■ (x + 16) 

(x + 1) • (x + 7) ■ (x + 17) 

(x + 2) ■ (x + 8)-(x + IS) 

. (x + 3) ■ (x + 9) ■ (x + 19) . 

fd T f 8 T 

#3: SOLVE — (x-(x + 6)-(x + 16)), x = x = - v x = -12 

^dx J 3 > 

' (x - 2)-(x - H)-(x - 26) ‘ 

(x - 2) ■ (x - 8) ■ (x - 18) 

(x - 2)-(x - S).(x - 10) 

3 

#4: VECTOR ([p(x , a, -2)], a, -3,3)= (x - 2) 

(x + 1) ■ (x - 2) ■ (x + 6) 

(x - 2)- (x + 4)- (x + 14) 

. (x - 2)-(x + 7)-(x + 22) . 


(d 




"H r- 

SOLVE — (( 
cdx 

; x - 2). 

.(X + 4) ■ ( 

x + 14)) 

i X = 

J \ 


#S : 


x 


v x = -10 
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Table V 


Family Function 

Roots 

Derivative 

Roots 

x 2 (x + 5)(x - 7) 

0, -5, 7 

2x(x-5)(2x + 7) 

0,5,-- 

2 

x 2 (x + 5)(x + 2) 

0, -5, -2, 

x(x + 5)(4x + 5) 

0,-4^ 

x (x + 5)(x + 9) 

o, -5, -9 

2x(x + 3)(2x + 15) 

0,-3,-— 

2 

x 2 (x + 7)(x + 10) 

0, -7, -10 

x(x + 4)(4x + 35) 

0,-4,-— 

4 

x 2 (x + 9)(x + 14) 

o, -9, -14 

x(x + 12)(4x + 21) 

21 

0,-12,-— 

4 

The Quintic 


(12) 

1 


Transformation 
(x ± k ) 2 (x + 5 a± k)(x - la + k ) 
(x ± k ) 2 (x + 5 a± k){x + 2 a + k) 
(x ± k ) 2 (x + 5a ± £)(x + 9a ± k) 
(x±k) 2 (x + l a± k)(x + 1 Oa ± £) 
(x ± k ) 2 (x + 9a ± £)(x + 14a ± £) 


In terms of multiple roots the quintic lends itself nicely to easy-to-work-with numbers which are 
small in quantity. However, for quintics of the form x(x-a)(x-Z?)(x-c)(x-<i), the derivative 
has no rational roots primarily because of Fermat’s Last Theorem whereby there are no integral values 
for which x 4 + y 4 = z 4 . Having run the computer through thousands of number combinations (just to 
be sure), no derivative with rational roots could be found. Our Table VI contains multiple roots only. 
Table VII contains the numerical families for the forms x 3 (x + a)(x + Z?) and x 2 (x + a) 2 (x + &). 


With all the patterns that do work, it was too tempting not to try to make a linear link among the 
cubic, quartic, and quintic forms. Let us examine the following facts: 


I. Cubic Form 

x(x - 5)(x - 8) 
Derivative 
(3x - 20)(x - 2) 

Smallest Root 
0 

2 

Largest Root 
8 

20 

3 

Range 

8 

14 

T 

Sum of Roots 
13 

26 _ 520 
3 ~ 60 

II. Quartic Form 

Smallest Root 

Largest Root 

Range 

Sum of Roots 

x(x - 3)(x - 4) (x - 7) 

0 

7 

7 

14 

Derivative 





2(x - l)(2x - 7) (x - 6) 

1 

6 

5 

21 _ 630 





2 ~ 60 


III. Quintic Form Smallest Root Largest Root Range Sum of Roots 
x(x - 2)(x - 3)(x - 4)(x - 6) 0 6 6 15 


Derivative 

5x 4 - 60x 3 + 240x 2 - 360x + 144 


0.616036... 


5.383963... 


4.767926 


12 
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I would never have attempted all this work without the user-friendly program DERIVE, as I was 
able to test many polynomials in seconds and quickly find derivative and corresponding factored 
forms. It goes without saying that the same Diophantine process can be applied to rational forms, 
making life a little easier for the curve sketcher. 


Part 2 of the Diophantine Polynomials will follow. 


Additional Editor’s Comments 

I enjoyed this contribution very much. As you might know from earlier articles I liked to use 
the computer as a training partner for the students providing random generated problems 
from various fields of the curriculum. 

Duncan mentioned that he used the TABLE function of the TI-83 for discovering certain ex- 
pressions as perfect squares. I felt inspired to use the Tl-Nspire for producing random gen- 
erated polynomial functions. Take the quartics presented in Table IV (the Pythagorean Trip- 
lets). 

First of all we can - preferable together with the students - derive the condition(s) for having 
rational roots of the first derivative. (These are the first three lines in the Calculator Applica- 
tion of NspireCAS below.) Then let’s discuss how to produce Pythagorean Triplets - and 
finally use the spread sheet application to create random Pythagorean Triplets. 

if x = m 2 -n 2 , y = 2 mn, z = m 2 + n 2 
then x 2 + y 2 — z 2 

This is shown in the next screen shot. 


g(x j a, b, c ) : =x • (x+ a ) • ( x+ b ) ■ (x+ a+ h ) 

Done 

■^-[g[x,a,b,c)) 

4'X 3 +6' 

(ft+Z?)'X^ + (2' £7^ + 6' Q'b+2' )'X+£7' (c7+Z?)' b 

zeros (g(x,, a , b, c ) )j xj 

\tQ^+b* —Q—b "(J a^+b^ +c?+£>] - (c7+£j) , 

to 

to 

to 

zeros f — (g(xj 6 0 } ~ 3 2j 2 8 ) )j x 
\dx 

{-48, -14, 20} 

zerosfygUlS, 1 2, 47) ),xj 


z er os(-^-(g(x, 8, -6, 2) ),xj 


{-6, -1,4} 


In columns A and B are defined random generated integer numbers between -10 and +10 
(without 0). Cell A1 and B1 contain the same formula which must be copied down - say 20 
times. 

In columns C, D and E are the values for the constants a, b, and a + b of the polynomial g. 


As you can see we obtain pretty much nice triplets of numbers. 
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A m 

Bn 

c a 

■b 

■ 


G 

1 

I 

F r 

♦ 

1 



=abs(a[] A 2-b[] A 2) 

=2*a[]*b[] 

=c[]+d[] 






8 

-2 

60 

-32 

28 






2 

-3 

-9 

72 

54 

126 






3 

1 

6 

35 

12 

47 






4 

-2 

-1 

3 

4 

7 






5 

6 

1 

35 

12 

47 






6 

-10 

8 

36 

-160 

-124 






7 

-3 

-4 

7 

24 

31 






8 

1 

5 

24 

10 

34 






9 

-1 

7 

48 

-14 

34 






10 

4 

7 

33 

56 

89 






11 

1 

7 

48 

14 

62 






12 

-7 

-6 

13 

84 

97 






13 

10 

3 

91 

60 

151 






14 

-4 

-9 

65 

72 

137 






15 

10 

-5 

75 

-100 

-25 






16 

-1 

3 

8 

-6 

2 






17 

7 

-10 

51 

-140 

-89 






18 

-7 

-5 

24 

70 

94 






19 

1 

7 

48 

14 

62 







CL 

Cl 

~7 ~7 

'fC 

J-±- 





5? 


A1 


= ra ndint(l,10)'(-l) randint ^- 2 ) 


< > 


The next three lines of the calculator present a test, if the roots of the first derivative are 
really rational numbers. I took the values of rows 1, 5, and 16. 

Let the students confirm the calculation by plotting one or more function graphs and investi- 
gating the turning points: 
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The last column of Duncan's Table IV gives advice how to create many polynomials from 
one “mother polynomial” applying easy to perform transformations. 

I add a small program (function) how to generate random quartics of this type having fixed 
the values a and b of the “mother function”. 


quartic(8,-6) produces “children of x (x + 8) (x- 6) (x + 2). 



Don’t forget defining a randseed when using the function another time. Otherwise you might 
find the same polynomials. You can extend this function by implementing a random gener- 
ated “mother”, too. 

Add more types of quartics and ... 

It should be no problem transferring this procedure to the Voyage 200 and/or the good old 
TI-92. 

This would not be the “DERIVE Newsletter J ’ if I wouldn’t add a DERIVE comment, too. Dun- 
can wrote in his paper that he used DERIVE to find many connections between the various 
types of polynomials and to quickly derive and factorize. 

Let me use DERIVE to create problems (again “members of certain families”) more auto- 
mated and to find representants of the families more systematically. 

I start with Table I, next to the last example. 

The first expressions confirm Duncan’s entries in his Table I 
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#16: S0LUTIQNS(cubl(7) = 0, x = - 


28 56 7 


5 5 10 


#17: SOLUTIONS — (250- x - 1575-x - 14700-x + 10976) = 0, x 

^.dx j 


400-x 2B6 


#1B: cubl - 


2B0-x + lBO-x - 


#19 : SOLUTIONS cubl 


0. x = 


1 & 16 


IB IB IB 


#20: SOLUTIONS — cubl - 

^.dx L 


o, X = - 


Luckily it worked. We can ask ourselves, how to achieve integer roots for the given function 
and the derivative as well? Inspecting #10 and #1 1 from above, we conclude that k must be 
a multiple of 10 (divided or not by 250). Ok, let me try: 


cuhlC-30) 


+ 27-x - 1080-x - 3456 


( cublC-30) 

#22: SOLUTIONS = 0, x = [-3, 24, -48] 

L 250 J 

fd cubl(-30) 'j 

#23: SOLUTIONS = 0, x = [12, -30] 

Ldx 250 J 

cubl(llO) 3 2 

#24: = x - 99 -x - 14520- x + 170368 

250 


#25: FACTO R(x - 99-x - 14520-x + 170368) = (x - 11) -(x + 88) ■ (x - 176) 

fd 3 2 'l 

#26: SOLUTIONS — (x - 99-x - 14520-x + 170368) = 0, x = [-44, 110] 

c.dx ^ 


Finally I am treating the family of family functions appearing in Table 


#27: — (x-(x + a)-(x + b)) = 3-x + x-(2-a + 2-b) + a-b 

dx 


#2S: S0LUTI0NSC3-X + x-(2-a + 2-b) + a-b = 0, x) 


-/(a - a-b+b) + a+ b -/(a - a-b + b) - a - b 


I asked myself how to find suitable “partners" which might be able to found a “family" to use 
also Duncan’s very charming notation. 

Then it would be necessary to have a function which presents possible partners to one given 
root, resulting in integer roots for the function and rational roots for the derivative. 
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#30: partner(a, 1, lj) := 


SELECT 


v = FLOOR (v ), v, VECTOR ([b, ^(a - a-b + b ) 
. 2 2 


], b, 1, u) 


Ul 


#31: partner(3 , -100, 100) = [-3, 0, 3, 8] 


The VECTOR-command produces pairs of b with I < b < u and the root - with a given a. The 
SELECT command filters all pairs with an integer root and collects only the first elements of 
the pairs - these are possible parters to found a family. 


#31: partner (3, -100, 100) = [-5, 0, 3, 8] 


#32: VECTOR 

SOLUTIONS 

fd 1 

— (x-(x + 3)-(x + p)) = 0, x 

, p, [-5, 0, 3, 8] 



dx j 

J 



#33: 


0 -2 
-1 -3 


-6 


4 

3 


There are only 4 possible partners between -100 and +100. (There are the same four 
between -10000 and + 10000 - after 55 seconds calculation time!). The roots of the derivate 
for all partners are given in #33. Let’s have another try. Possibly 7 would be more attractive 
for partners? 

#38: partner(7, -100, 100) = [-33, -8, 0, 7, 13, 40] 


#39 : VECTOR 

SOLUTIONS 

fd ^ 

— (x-(x + 7)-(x + p)) = 0, x 

, p, [-33, -8, 0, 7, 15, 40] 


V 

^.dx > 

J 


21 


11 

3 


-4 


14 

3 


#40: 


0 


14 

3 


-7 


7 


3 


-3 


33 

3 


-28 


10 

3 


My Conclusion: I felt very inspired by Duncan’s contribution. I find it useful mainly because 
of two reasons: it makes math teacher’s daily work easier and it presents a lot of problems 
which could (and should) discussed with the students, e.g. filling in the Tables’ last columns 
which are “Conditions for integral or rational roots” and the possible “Transformations” or 
including considerations about the 2 nd derivative, ... And this is much more than only a tool 
for generating bundles of “beautiful” investigations of polynomial functions, Josef. 

See also Duncan’s mail on page 38. 
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Ein Mathematisches Modell fur Schneckenhauser 

A Mathematical Model for Snail Shells 

von Piotr Trebisz 

Piotr-Andrze j .Trebisz6gmx.de 

Wahrscheinlich ist schon jedem DERIVE -Nutzer der spiralformige Toms aufgefallen der auf der Ver- 
packung von Derive abgebildet ist. Dieses geometrische Objekt hat mich dazu inspiriert, einen Beitrag 
uber Schneckenhauser zu schreiben, die auf logarithmischen Spiralen basieren. 

Schneckenhauser sind ein sehr schones Beispiel dafur, dass logarithmische Spiralen ein haufig wie- 
derkehrendes Muster in der Natur sind. Logarithmische Spiralen sind selbstahnlich, das heifit, dass 
eine kleine Spirale genauso aussieht wie eine grobe Spirale. Das ermoglicht Schnecken ihre Hauser 
kontinuierlich weiterzubauen, ohne dass sich deren Geometrie andert. 

Ich werde im Laufe der folgenden Artikelreihe ein mathematisches Modell fur solche Schneckenhau- 
ser vorstellen. Da die von mir vorgestellten Schneckenhauser einige besondere Eigenschaften haben, 
die so in der Literatur noch nirgendwo beschrieben wurden, nehme ich mir die Freiheit heraus, sie 
nach mir zu benennen: Die Trebisz-Spiralen. 

Beginnen wir bei der Konstmktion unserer Schneckenhauser mit einer einfachen Variante, bei der uns 
als Basis eine planare logarithmische Spirale dient, ahnlich wie bei einer Nautilus-Schale. Ich nenne 
sie „Trebisz-Spirale - Typ 1 “ 



It is very likely that all DERIVE users have noticed the spiral formed torus on the package 
and the manual. This geometric object inspired me to write a paper on logarithmic based 
snail shells. 

Snail shells are a very beautiful example for the fact that logarithmic spirals are a frequently 
recurring pattern in nature. Logarithmic spirals are self similar, i.e. a small spiral looks like a 
large one. And this enables snails building their shells continuously without changing their 
geometry. 
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Over the following series of contributions I will introduce a mathematic model for such snail 
shells. As these snail shells will show some special properties, which I didn’t find elsewhere 
in literature, I take the freedom to call them Trebisz Spirals. 


We will start constructing our snail shells with an easy model, taking a plane logarithmic spi- 
ral as base, very similar to a shell of a nautilus. I call it “Trebisz-Spiral - Type 1”. 


Als Geriist fur das Schneckenhaus dient eine planare logarithmische Spirale die nach ihrem Radius 
bzw. ihrer Lange parametrisiert ist. Der Parameter .,/" steht hierbei fur den Radius bzw. die Lange der 
Spirale, wahrend ..h" die Basis des Logarithmus ist. 

A plane logarithmic spiral which is parameterized according its radius and its length, respec- 
tively is serving as the shell’s scaffold. Parameter T is standing for the radius and the length 
of the spiral, respectively, and “£>“ is the base of the logarithm. 


#1: [CaseMode := Sensitive, InputMode := Word] 

#2: [b :e Real (1, «) , m :e Real [0, «) , t :e Real [0, «) , s :e Real [-n, +rr]] 

#3: [w :e Real, n :e Integer (0, °o) , k :e Integer] 


#4: SNAIL (b, t) := [0, t-C0SC2-rr- L0G(t, b)), t- SIN(2-rr- L0G(t, b))] 


Exampl e : 

SNAIL(e, t) SUB [2,3] 
wi th 

0 < t < 1 



The frequent reader of the DNL knows 
my devotion to sliders. We can ani- 
mate the base spirals in the 3D-plot 
window: 

SNAIL (bj 


See the respective plot parameters to 
for presenting 5 rotations: 



Minimum 

Maximum 

Number of Panels 

s: 

|£ 

b 

b 

t: 

1° 

l 5 

|200 



Josef 
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Dass diese Spirale nach ihrem Radius parametrisiert ist, lasst sich leicht uberpriifen, indem man sich 
anschaut, wie weit ein Punkt auf dieser Spirale vom Zentrum entfemt ist in Abhangigkeit von 

It is easy to check that the radius is the parameter of this spiral by inspecting the distance of 
a point on the spiral to the centre dependent on T. 

#S: | SNAIL (b, t) | = t 

Dass die Spirale auch nach Ihrer Lange parametrisiert ist, lasst sich dadurch zeigen, dass die Tangente 
an jedem Punkt der Spirale eine konstante Lange hat 

The fact that this spiral has also its length as parameter can be shown by demonstrating that 
the tangent on each spiral point has constant length. 

2 2 

d -/(LN(b) + 4 -tt ) 

#6: — SNAIL (b, t) = 

dt LN(b) 

Wie man sieht, ist die Spirale auch nach der Lange parametrisiert, aber multipliziert mit einem kon- 
stanten MaBstabs-Faktor. 

Um den Mantel des Schneckenhauses zu konstruieren benotigen wir ein lokales mitdrehendes Koordi- 
natensystem aus 2 Vektoren die senkrecht zueinander und zur Tangente der Spirale stehen. 

Dies gelingt uns, indem wir die 2. Ableitung der Spirale nach ,J" bestimmen. Dabei machen wir uns 
zu Nutze, dass die 2. Ableitung einer Kurve, die nach ihrer Lange parametrisiert ist, einen Vektor er- 
gibt, der in der Kriimmungsebene der Kurve liegt und der senkrecht zu ihrer Tangente steht. 

As one can see the spiral is also parameterized wrt its length multiplied by a constant scaling 
factor. 

For constructing the surface of the snail shell we need a local co-rotating system of coordi- 
nates consisting of two vectors which are perpendicular to each other and perpendicular wrt 
the tangent of the spiral. 

The 2 nd derivative of the spiral wrt T can help. With its length as parameter the 2 nd derivative 
of the curve results in a vector lying in the osculation plane of the curve perpendicular to its 
tangent. 

d ^ fd ^2 

#7: — SNAIL (b, t) . — SNAIL (b, t) = 0 

cdt J Ldt; 

Das Kreuzprodukt aus Tangente und 2. Ableitung liefert uns einen weiteren Vektor der senkrecht zur 
Tangente und zur 2. Ableitung steht. Dieser Vektor soil unsere lokale mitdrehende X-Achse sein. Und 
die lokale Y-Achse konnen wir direkt aus der 2. Ableitung nach „t“ bestimmen. Die lokale mitdrehen- 
de Z-Achse entspricht der Tangente der Spirale. 

The cross product from tangent and 2nd derivative delivers another vector which is perpen- 
dicular to both tangent and 2nd derivative. This vector shall form our local co-rotating X-axis. 
The local Y-axis is given directly by the 2nd derivative wrt “t”. the local co-rotating Z-axis is 
the spiral tangent. 
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Yd 1 


fd 1 

2 ") 

#S : X_Achse(b, t) := SIGN 

— SNAIL (b, t) 

X 

— 

SNAIL (b, t) 

\ 

.^dt j 


^dt; 

J 



Yd " 

2 

#9: Y_Achse(b, t) := - SIGN 

— 

SNAIL (b, t) 


-:dt_. 

j 


# 10 : 


Z_Achse := SIGN 


fd 

— SNAIL (b, 
jJt 


t) 


Es lasst sich zeigen dass dieses lokale mitdrehende Koordinatensystem die selbe Handigkeit besitzt 
wie das globale Koordinatensystem, in dem die Spirale defmiert ist. 

We can show that this local co-rotating system of coordinates has the same direction as the 
global system where the spiral is defined. 

#11: X_Achse(b, t) x Y_Achse(b, t) = Z^Achse 


Die Mantelflache soil einfachheitshalber einen kreisformigen Querschnitt haben. Der Radius der Man- 
telflache soil linear mit dem Radius bzw. mit der Lange der Spirale wachsen, „m“ ist der Linearfaktor. 

To make it easy we will have a circle as cross section of the surface. The radius of the sur- 
face shall increase linearly according to the radius - and the length - of the spiral with "m" as 
linear-factor. 

#12: Mantel (b, rn , t, s) := m- 1- (COS(s) ■X_Achse(b l t) + SIN(s) ■ Y_Achse(b, t)) 


Und schon ist das Schneckenhaus fertig. Es ist die Summe aus Spirale und Mantel. 

Here we go, the snail shell is ready. It is the sum of spiral and surface. 


#13: Spi ra!TorusYZ_R(b, m , t, s') Mantel (b, m , t, s) + SNAIL (b, t) 


#14: Spiral To rusYZ_R 


0.5, , t, s 

3 y 


0 < s < 271, 0 < t < 1 ; 40 panels 
each 



Will man das Schneckenhaus nicht nach dem Radius, sondem nach Zahl der Umdrehungen paramet- 
risieren, dann ersetzt man „t“ durch ..h i". Der Radius wachst dann nicht mehr linear zu sondem zu 
b‘. 

When we don’t take the radius as parameter but the number of rotations then we have to 
substitute T by "b A t". Then the radius does not increase linear wrt T but to b\ 
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#1S: 


t 

SNAIL (b, b ) 


t 


= b 


t t 

#16: Spi' ralTorusYZJCb, m , t, s) := Mantel (b, m, b , s) + SNAIL (b, b ) 

Here again the slider bars give some useful insight about the importance of the parameters b 
and m. 

I introduce the sliders in order to animate the “linear-factor m” and the base of the spiral as 
well. 

See the snail shell with t = 2 and two combinations of b and m. 


Josef 



0<t<l 


-1 <*< 1,2 


Es gibt leider ein Problem mit dem Li- 
near-Faktor „m“. Wird er zu groB ge- 
wahlt, dann durchdringt sich die Man- 
telflache bei jeder Umdrehung selbst. 
Beispiel: b = e und m = 2/3. 


#17: 


Spi' ralTorus'YZ_R 


e. 



0 < t < 1 and -n < s < n 


Unfortunately we are facing a prob- 
lem with the linear-factor "m". 
Choosing m too big then the surface 
will interpenetrate at each rotation, 
e.g. b = e and m = 2/3. 
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Wird „m“ zu klein gewahlt, dann klafft 
bei jeder Umdrehung eine hassliche 
Liicke zwischen der Mantelflache. Bei- 
spiel b = e und m= 1/4. 


# 1 £: 


Spi' ra1TorusYZ_R 


1 ' 
4 , 


Choosing "m" too small, then we 
have an ugly gap between the sur- 
face, e.g. b = e and m = 1/4 



Wie kann man zur Basis „b“ das richtige ..m" so bestimmen dass sich die Mantelflache nach jeder 
Umdrehung exakt beriihrt, ohne sich selbst zu durchdringen? b = e und m = 0.4593761774. 

How to determine the right “m“ for given base “b“ that the surface is osculating exactly after 
each full rotation without interpenetrating? b = e and m = 0.4593761774. 



Dazu mussen wir uns den Schnitt des Schneckenhauses durch die XY-Ebene anschauen. Die Z-Koor- 
dinate von #16 (vereinfacht) muss also 0 sein. Wir losen die entsprechende Gleichung nach „t“ auf: 

We have to investigate the section of the snail shell with the XY-plane. Hence the 
Z-coordinate of simplified #16 must equal zero. We solve the respective equation for T: 


# 21 : 


2-TT-m- SIN(s) 


+ 1 


y(LN(b) 


2 

4-tt ) 


■ SIN(2-T7-t) 


m- LN(b) ■ SIN(s) -COSC^-tt- t) 


2 2 
^/(LN(b) + 4-rr ) 


= 0 


f m- LN(b) ■ SIN(s) > 

ATAN 

2 2 

#22: l yCLN(b) + 4 -tt ) + 2-rr-m- SIN(s) J 

t = 

2-tt 

Wir konnen die gefundene Losung verallgemeinern, denn das Schneckenhaus schneidet die XY-Ebene 
nicht nur einmal, sondern bei jeder halben Umdrehung. Wenn wir also zu der Losung ein beliebiges 
ganzzahlige Vielfaches von 1/2 dazu addieren, dann erhalten wir wieder eine giiltige Losung. „k“ ist 
eine beliebige ganze Zahl. 


We can generalize the found solution because the snail shell does not intersect the XY- 
plane only once but at every half rotation. We obtain another valid solution by adding any 
integer multiple of 14. ”/c” is an arbitrary integer. 
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ATAN 


rn- LN(b) ■ SIN(s) 


# 22 : 


2 2 

-/(LN(b) + 4-n ) + 2-TT-m- SIN(s) j 


2-n 


Fur die Losung des Selbstdurchdringungsprolems konnen wir aber nur jene Losungen fur gebrau- 
chen, die sich in einer ganzen Umdrehung unterscheiden, deshalb schranken wir die Losungen fur „t“ 
wie folgt ein: 

For solving the interpenetrating problem we can use only such solutions for T which differ in 
a full rotation. We set a constraint for „f as follows: 


#23: 


t = 


ATAN 


2 

l J (LN(b) 


m- LN(b) ■ SIN(s) 


2 

+ 4 -n ) + 2-TT-m- SIN(s) j 


2-tt 


+ k 


Setzt man diese Losung fur in die Schneckenhaus-Formel ein, wird die Z-Koordinate 0 und man 
erhalt den XY-Schnitt der Flache, wo man sehen kann, ob sie sich selbst durchdringt oder nicht. 
Schauen wir uns mal ein paar Schnittkurven fur k = 0 an: 

Substituting this solution for T in in the formula for the snail house the Z-coordinate will be- 
come 0 and one obtains the XY-intersection of the surface. This section shows whether the 
surface is interpenetrating or not. 

We will look at some intersection curves for k = 0: 




m- LN(b) ■ SIN(s) ^ 




ATAN 






2 2 





< ^(LN(b) + 4-n ) + 2-TT-m- SIN(s) j 



SUB ST 

Spi' ra1Toru5VZ_U(b, m, t, s'), t, 




X 


2-n 

J 

J 


[ 1 , 2 ] 


#2B: 


2 

ATAN (m ■ LN (b) ■ SIN (s) / U (LN (b) 

m- b 


2 

+ 4 -n ) + 2-TT-m- SIN(s)))/(2-tt) 

■cos( s ), 


2 2 

ATAN (m ■ LN (b) ■ SIN ( 5 ) / (./(LN (b) + 4-tt ) + 

b 


2 2 2 2 2 2 
2-7r-rn-SIN(s)))/(2.7r) SIGN(^(LN(b) + 4-n ) + 2-TT-m- SIN(s)) -^(^(LN(b5 + 4-n ) ■ (m -SIN(s) + 1) + 4-7r-rn-SIN(s)) 


2 2 1/4 

(LN(b) + 4-n ) 


b = e and m = 3/2 


-l 


l 



b = e and m = 1 
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Hoppla!!! Warum ist denn im 2. Beispiel die Schnittkurve unterbrochen? Der Grund dafiir liegt in der 
Losung, die wir fiir „t“ bestimmt haben. Formal gesehen ist die Losung vollkommen in Ordnung. 
Wenn man sie zur Probe in die Z-Koordinate des Schneckenhauses einsetzt, dann wird diese 0. 

Nur hat die Losung einen gravierenden Schonheitsfehler! Bei bestimmten Kombinationen von „Z>“ und 
„/«“ zerfallt die Schnittkurve in zwei Teile, die zwar genau zusammenpassen, sich aber auch in genau 
einer halben Umdrehung unterscheiden! Das erkennt man daran, dass sich das zweite Stuck der 
Schnittkurve auf der anderen Seite der X-Achse befindet. Dieses Kurvenstiick gehort zum Schnitt des 
Schneckenhauses mit der XY-Ebene nach der nachsten halben Umdrehung! 

Offensichtlich gibt es einen bestimmen Wertebereich von ,.v", wo es in der Y-Koordinate zu einem 
Vorzeichenwechsel kommt. Und das wiederum fiihrt zu einem „Umdrehungswechsel“ in der Schnitt- 
kurve! Dieser „Umdrehungswechsel“ zeigt sich sehr schon, wenn man sich die Schnittkurven fur zwei 
Losungen anschaut, die sich in einer halben Umdrehung unterscheiden. 

Oops!! Why is there an interruption in the 2 nd example? The reason lies in our solution for T. 
Formally the solution is absolutely ok. The Z-coordinate is becoming 0 by substituting this 
solution, but ... 

... it is a solution with a serious flaw. Certain combinations of “b“ and “m“ result in a decom- 
posing of the intersection curve in two parts which - fitting exactly to each other - differ in 
exact one half rotation! This can be recognized by the fact that the second part of the inter- 
section curve is lying on the other side of the X-axis. This part of the curve is part of the in- 
tersection of the shell with the XY-plane after the next half rotation! 

Obviously there is a certain domain of "s" which results in a sign change of the Y-coordinate. 
And this leads subsequently to a change of the rotation in the intersection curve! This “rota- 
tion change” appears very clear by observing the intersection curves for two solutions which 
differ in a half rotation. 


ti 


ATAN ■ 


2 

L yCLNCb) 


m ■ LN(b) ■ 5IN(s) 


2 

+ 4-r r ) + 2 ■ Jr ■ m ■ SIN(s) ^ 


2 -7T 


und/and 


t2 


f 

ATAN ■ 


2 

L 7(LN(b) 


m ■ LN(b) ■ SIN(s) 


2 

+ 4-7T ) + 2 ■rr-m ■ 5IN(s) > 


2 -7T 


1 

2 


Wie man sieht, fiigen sich die Teile der zwei 
Schnittkurven passgenau zusammen. Gibt es eine 
Moglichkeit die Losung zu reparieren? 

Es gibt tatsachlich einen Weg, dazu mussen wir 
uns genauer anschauen fiir welche Wertebereiche 
von „s“ der Vorzeichenwechsel stattfmdet. 
Gucken wir uns dazu die 4 Terme genauer an, 
aus denen die Y-Koordinate der Schnittkurve 
zusammengesetzt ist. 


b = e and m = 3/2 
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Nachdem es uns endlich gelungen ist die Schnittkurve des Schneckenhauses mit der XY -Ebene pro 
Umdrehung zu berechnen, konnen wir uns daran machen das Selbstdurchdringungsproblem zu losen. 
Dazu miissen wir den Parameter „m“ so bestimmen dass sich die Schnittkurven von zwei aufeinander 
folgenden Umdrehungen exakt beriihren ohne sich zu schneiden. Schauen wir uns die Schnitte fur 
k = 0 und k= 1 an. 

Substituiere fxir „t“ aus #33 (fur k - 0 und dann fur k = 1) in SpiralTorusYZ U und setze dann die 
entsprechenden Wertepaare fur „Z>“ und „m“ ein. (Substituiere wie in #24 um #34 und #35 zu erhalten 
- hier aus Platzgriinden nicht angegeben.) Das ergibt die folgenden Schnitte. 

Having succeeded calculating the intersection curve of the shell with the XY-plane per rota- 
tion we can start solving the interpenetration problem. We have to determine the parameter 
“m“ in such a way that the intersection curves of two consecutive rotations are osculating 
without intersecting. Let’s inspect the intersections for k = 0 and k = 1 . 

Substitute for T from #33 (with k = 0 and then with k= 1) in SpiralTorusYZ_U and then sub- 
stitute again the respective values for b and m. (Substitute like in #24 to obtain #34 and #35, 
which are both omitted here to save space.) This gives the following intersections: 

b = e,m = V 4 b = e, m = 312 b = e, m = 0.4593761774 



Wie man sieht durchdringen sich die beiden Schnitte entlang der Y-Achse auf gegenuber liegenden 
Seiten der Schnittkurven. Die X-Koordinate muss also 0 sein. Wir losen nach „s“. 

As you can see the both intersections interpenetrate along the Y-axis at opposite sides of the 
curves. Hence, the X-coordinate must become 0. We solve for "s". 
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#37: SOLUTIONS Ui-b 


2 2 2 
(SIGN(,/(LN(b) + 4-71 ) + 2.ir m-SIN(s)) - l)/4 + ATAN(m- LN(b) -SIN(s)/ (7(LN(b) + 


4-77 ) + 2 -77-m-SIN(s)))/(2 -77) r 77 77 3 -77 

• COS(s) :0,S J = , , 

.2 2 2 

(22 2 

(SIGN(7(LN(b) + 4-77 ) + 2 *77 • m • SIN(s)) + 3)/4 + ATAN(m- LN(b) -SIN(s)/(#(LN(b) + 


#38: SOLUTIONS U-b 


4-77 ) + 2 -77- 111 • SIN (s)))/ (2 -77) I 

77 

77 

3-77 

-COS(s) = 0 , sj = 

, 

— i 



2 

2 

2 


Die Losungen fur a sind: / The solutions for s are: 

71 77 

k = 0 A s = and k = 1 a s = - — 

2 2 

Die beiden Losungen werden eingesetzt (in #35 und #36). 

Both solutions are substituted (in #35 and #36). 


#39: 0, b 


#40 : 0 , b 


2 2 2 2 2 
ATAN(m- LN(b)/(#(LN(b) + 4-77 ) + 2 •77*ra))/(2 *77) #((m + l)-</(LN(b) + 4-77 ) + 4-77-m) 

b 

2 2 1/4 

(LN(b) + 4-77 ) 

2 2 2 2 
(3 + SIGN(#(LN(b) + 4-77 ) - 2 • 77 - m))/4 - ATANfm • LN(b)/(#(LN(b) + 4-77 ) - 


2 2 2 

2 -77 -in))/ (2 *77) </((m + l)-#(LN(b) + 4-77 ) - 4-77-m) 

2 2 1/4 

(LN(b) + 4-77 ) 

Wenn sich die Mantelflache nur beriihren, aber nicht durchdringen soil, miissen beide Y-Koordinaten 
gleich sein. Folgende Gleichung muss also nach „m“ gelost werden: 

zweite Komponenten von #39 = zweite. Komponente von #40. 

If the surface shall only osculate (and not interpenetrate) then both Y-coordinates must be 
equal. So we have to solve the following equation for “m“: 

second component of #39 = second component of #40. 


#41: b 


2 2 2 2 
(3 + SIGN(#(LN(b) + 4-77 ) - 2 -77-m))/4 - ATAN(m - LN(b)/(#(LN(b) + 4-77 ) - 


2 2 2 2 2 

2 -77 • m))/(2 -77) 7((m + l)-#(LN(b) + 4-77 ) - 4-77-m) ATAN(m- LN(b)/(#(LN(b) + 4-77 ) + 

= b 

2 2 1/4 

(LN(b) + 4-77 ) 

2 2 2 

2 -77 -m))/ (2 -77) >/((m + l)-7(LN(b) + 4-77 ) + 4-77-m) 

2 2 1/4 

( LN (b ) + 4-77 ) 
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Gesucht ist ein Wert fur „m“ fur den gilt t'(h.m ) = 0. Fiir das Newton-Verfahren wird die 1. Ableitung 
von f '(b,m) nach „/h“ benotigt. g (b,m) ist die 1 . Ableitung von f {b,m). 

We would like to find a value for "m" with f (b,m) = 0. For the Newton-algorithm we need the 
first derivative of f (b,m) wrt "m". g(b,m) is the first derivative of i(b,m). 

d 

#46: g (b , m) := — f(b, m) 

dm 

Leider ist DERIVE nicht selbststandig in der Lage die 1 . Ableitung von f (b,m) nach „/n“ zu bestim- 
men. Jeder Versuch fuhrt zu einer Fehlermeldung wegen zu wenig Speicher. 

Wenn wir aber DERIVE exakt vorgeben wie es die 1 . Ableitung bestimmen soil, ist die Software doch 
dazu in der Lage. Dazu zerlegen wir i(h.m ) in 4 Teilfunktionen . . . 

Unfortunately DERIVE is unable to find this first derivative. Each attempt leads to an error 
message: Memory Exhausted. 

If we give some hints how to find the first derivative then the software is able to differentiate 
f (b,m). We decompose f (b,m) into 4 partial functions ... 


(I omit printing this - from the didactical point of view - interesting part of Piotr 's file in order 
to save space. I invite you to follow his steps in expressions #47 - #55 of his dfw-file. #55 is 
the final result for g (m,b). See also my comment at the end of this article, Josef) 

2 2 2 2 
(3 + SIGN(V(LN(b) + 4-7? ) - 2-71 •iii))/2 - ATAN(lH- LN(b)/G/(LN(b) + 4-ff ) - 

#55: g(b, m) := b 


2 •«))/>? 


2 2 LN(b) 

2-m-V(LN(b) + 4 -n ) - - 4-ji 


2 2 


- b 


ATAN(m-LN(b )/(■/( LN(b) + 4-/7 ) 


2-jt-ii))/ji 


2 2 LN(b) 

2-m-/(LN(b) + 4 -n ) + + 4-77 


Den Newton-Algorithmus kann man als kleines DERIVE - Programm realisieren. 

The Newton-Raphson algorithm is realised in form of a short DERIVE program. 


h(b, ni, ii ) := 

Prog 

6m :g Real 

#56: WRITECAPPENDC'm : = " , APPROX(STRING(m) , n))) 

6m := m - APPROXCfCb, m) , 2 .n)/APPROX(g(b , in), 2-n) 
RETURN APPROX (6m, n) 


M_YZ(b , ii ) := 

Prog 

#57: V :e Vector 

V := ITERATES (h(b , in, n), m, 0) 
RETURN ViDIH(V) 
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M_YZ(b, n) berechnet das passende ”m" zur Basis "b" auf "n" Stellen genau. Und da wir auch die 
Windungszahl berechnet haben, konnen wir nun das Schneckenhaus auch nach der Windungszahl 
parametrisieren. Das findet sich in SpiralFlaecheYZW. 

M_YZ(b, n) calculates the matching "m" fort he base "b" accurate to "n" decimal places. And 
as we have calculated the coil number, too, we are able to parameterize the snail house ac- 
cording to this number, giving SpiralFlaecheYZ_W. 

Now we can find the answer for my question on page 23: 

“How to obtain this very special m?” 

2071448 

#63: H_YZ(e, 10) = 

4509263 

#64: H_YZ(e, 10) = 0.4593761774 

Here I came across a strange - bug? Preparing the publication of this article I reproduced all 
calculation steps and I failed simplifying (and approximating as well) expression #63. I asked 
Piotr and he answered that he is aware of this fact, but there will occur nor problems in the 
respective mth-file. So I stored his dfw-file as an mth-file - loosing all graphs and comments - 
and surprisingly enough, now it worked - as you can see above. 

What happened here? I switched back to the dfw-file and opened the function g(m,b) via 
Author > Function Definition. I saw the function for a moment and then there was in place of 
the function the message “Memory Full” in the Edit Window!! I tried calculating the first 
derivative of f(m,b) - no success, Memory Exhausted. 

When I did the same in the mth-file, the same message appeared in the Function Definition 
Window, but when simplifying M_YZ which is based on f(m,b) and the predefined g(m,b) 
DERIVE continues cooperating with the user ... These are expressions #63 and #64 above. 
Don’t ask my why. Hopefully Albert Rich is reading this remark and maybe that he will 
provide some explanation for this strange behaviour. 

If one of our members knows the reason or if he has some conjectures then please let us 
know, Josef (and Piotr, of course) 

I continue with Piotr’s closing lines. 

Ich bevorzuge es immer auf n + 2 Stellen genau zu rechnen. Eine Extra Vor-Komma-Stelle und eine 
Extra Nach-Komma-Stelle fiir den Rundungsfehler 


I prefer calculating for n + 2 decimal places. One extra digit in front of the decimal point and 
one extra digit at the end of the number to consider a possible rounding error. 


#58: 

Sp"iralFlaecheYZ_R(b, n, 

t, s) := Spi ralTorusYZ_R(b , 

M_YZ(b , 

n) , 

t, s) 

#59: 

SpHralFlaecheYZJJCb, n, 

t, s) := SpH ralTorusYZ_U(b , 

M_YZ(b , 

n), 

t, s) 

#60: 

Sp"iralFlaecheYZ_W(b, n, 

t, s) := Sp"i ralTorusYZ_W(b , 

M_YZ(b , 

n), 

t, s) 
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Wir schlieBen mit ein paar hiibschen Schneckenhauser fur die Basen 2 und 6. 

We close presenting some nice snail shells for bases 2 and 6. 

Base 2 Base 6 



In Snail_shells.dfw finden sich alle Details. Die Zusammenfassung aller wichtigen Funktionen fur das 
Plotten der Schneckenhauser und fur weiteres Experimentieren bildet die mth-Datei 
Logarithmic_Spiral_Surfaces.mth. 

The respective dfw-file is called Snail_shells.dfw. All important functions for plotting the 
shells and experimenting with various parameters are collected in an mth-file, called 
Logarithmic_Spiral_Surfaces.mth. 

If there are any serious mistakes in the English translation - I am quite sure that there are 
some - then please don't blame Piotr Trebisz. I received the German version and tried to 
produce the English translation. 


Josef 
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As you can see in the following screen shots, MuPad and MATHEMATICA can both derive 
function f(b,m). The derivative is much bulkier than Piotr’s result. I plotted f (e,m) and g(e,m) 
in MuPad and compared with the DERIVE plot. It looks pretty the same, so the results should 
be identical? 




MuPad above and MATHEMATICA below 


ln[3]:= D [f [b, m] , m] 




m bg(6) 


_ 


0ut[3] : 


|= -(^2 'sj log 2 (i)+4^ 2 m 4- 4 ;rj 


Wnm+^j bgp(^)-t4^ j 


log(b) (4nm + (m 2 +1) V log 2 (&) + 4 n 2 ^ 


log <30 


2mxhg(b) L 

^2jirmW log^-H** 2 [2 jrm+V log 2 {b)+An 2 ^ ) 


tan-1 

1 htfm+Vlog^)-^* 2 
\h * 


m 2 log 2 (6) 


- + 1 


^ logoi' )- h4^ ) 

-if mUt(b) 


\og(b) | '(m 2 + 1) -\f \o^{b)-\-An 2 - A mii^ 


LN (b) 


2wigLN|ft) 


V log 2 {b)+A?t 2 -2 mn (V log 2 {b)+An 2 -2 m n} 


m 2 Ul(b) 2 


+ 1 


^ (V log 2 (30 +4 jt 2 -2 m jtJ j 


- sgn'^i/ log 2 ( b ) + 4 k 2 - 2 n- CenterDot ( ' 0,1 \/r, m) 


^-iL-aa — 

^ j (sgn(V bg 2 (^)44 a 2 -2 am) +3) log^ftMa 2 2 mn 
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Exchange of mails connected with Duncan McDouqal 's "Diophantine Polynomials" 


Mail from Duncan McDougal as answer on my request concerning the table in (12) on page 13. Dun- 
can had 740/60 as sum of the roots in the last row . . . 

In an earlier mail he had asked: “Just out of curiosity (if I may ask), what are your intentions with this 
article?” 


This was my answer: 

I have two intentions: 

(1) Finding and exploring patterns is very important in mathematics education and in 
mathematics in general. Your paper is an excellent example for discovering - 
and using - patterns. 

(2) The second reason is given in my final comments (page 21). 

Hi Josef, 

You are correct, the sum is 12 and should remain 12. However, the pattern becomes 
520, 630, and 740 as an arithmetic sequence and 740/60 reduces to 37/3, a wee bit more 
than 12. 1 am still secretly hoping that someone sees a pattern I don't see and will come up 
with a brilliant discovery on the Quintic. Thus, I am suggesting that although 12 is the prob- 
able answer, maybe 37/3 is what it could be??? I acknowledge that I've proven that there is 
no monic quintic for which the first derivative contains integral and/or rational roots, but the 
pattern is hard to resist or ignore. 

I am still very much interested in this topic and will entertain questions or comments 
from any interested party. I've only just begun the process for non-monic polynomials and 
already the quartics look daunting. 

Anyway, I am flattered that you like my work enough to use it because that is what it 
was meant for. 

I remain 

Yours truly, Duncan McDougall 


Danny Ross Lunsford (22 01 2011) 

Hello Derivers, the routines provided in the LinearAlgebra.mth file leave something to be desired in 
practice - has anyone resolutely attacked the problem of diagonalization and determination of eigen- 
spaces of large matrices? For my application large = 8x8. I need a robust and fast diagonalization 
scheme. (DNL#80). 

Hi Derivers, (19 02 2011) 

In implementing the Johnson Trotter algorithm one comes to a point where the direction of a mobile 
integer must be changed. My working area is an (N+2) x 2 matrix with the working permutations in 
the 1st column and the directions in the 2nd (directions are + or - 1). Eventually one has to change 
directions for all integers larger than the largest mobile integer - so I have a statement to run over the 
array 

MAP(IF(P SUB i SUB 1 > Largest, P SUB i SUB 2 := -P SUB i SUB 2), i, [2,..,N+1]) 
However the sub-array assignment never takes place, even though the statement is executed. What 
exactly is going on? I know sub-array assignments must be simplified, but that should be happening 
here. Thanks in advance! 

-drl 
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Hi Derivers, (19 02 2011) 

OK see attached - 1 changed the map statement a bit but still no joy. The SHAPE function is one I 
made that creates arrays, don't worry about it, the size of the array created is (N+2) x 2. The output is 
included and you can see that the sign changes are never executed. 

(20 02 2011) 

It doesn't work - see attached - a function to generate all permutations of 1 to N by accretion was given 
in Derive Newsletter 41. It creates two nested loops, the inner one with k going from m to 1 in steps of 
-1, the outer with m going from 2 to N in steps of 1. All the variables are local to the PERMS function. 
We can make the inner loop into a MAP statement but the syntactically identical function works when 
the accretion variable t is global, but fails when you make it local to the function - I'm completely 
baffled 

#1: Generate all permutations of N integers by accretion - from DNL41 


PERMSfN, k, m := 2, s := [[]], t) := 

Loop 

If m > N 

RETURN SORT (s) 
k := m 
t := [] 

#2: Loop 

If k = 0 exit 

t := APPENDft, VECTOR (I NSERT(m, v, k) , v, s)) 
k := -1 
s := t 
m :+ 1 


TPERMS(N , k, m := 2, s := [[1]]) := 

Loop 

If m > N 

RETURN SORT (s) 

#3: t := [] 

MAP(t := APPENDft, VECTOR (I NSERT (m , v, k) , v, s)), k, m, 1, -1) 
s := t 
m :+ 1 


#4: t := [] 


#5 : We r ve replaced the inner loop over decreasing k by a MAP statement 

#6: t is a global variable which we set to empty 

"12 3 

13 2 

2 13 

#7: TPERMS(3) = 

2 3 1 

3 12 

_ 3 2 1 

#S: It works - now make t a local variable to TPERMS 

TPERMS_(N , k, m := 2, s := [[1]], t) := 

Loop 

If m > N 

RETURN SORT (s) 

#9: t := [] 

MAPCt := APPENDft, VECTOR(INSERT (m , v, k) , v, s)), k, m, 1, -1) 
s := t 
m :+ 1 

#10: TPERMS_(3) = [[3, 2, 1]] 

#11: Doesn't work any more - why?? 
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Danny Ross promised an article on Clifford Algebra and Dirac matrices some time ago. I 
asked him whether his request is in connection with this intended paper for the DNL and he 
answered. 


20 02 2011 

Yes, will work out solution to hydrogen fine structure as an example application. See attached for 
teaser on the algebraic environment. Derive makes it very easy to keep track of the rather complex 
algebraic operations in the Clifford algebra of Dirac matrices. 


On page 48 in DNL#80 I published the email communication with Albert Rich. Unfortunately I 
did not include the respective website where you can find many details about his fascinating 
current research. Albert wrote: 

Hello Josef, 

Yes, feel free to include my 2 January 2011 email in the next DNL. Also please include the 
website address 

http://www.apmaths.uwo.ca/RuleBasedMathematics 

As far as my current research, I recommend the following ongoing discussion on the Usenet 
group sci. math. symbolic 

http://qroups.qooqle.com/qroup/sci.math.symbolic/browse thread/thread/90b4f2e39ee5ad9d 
?hl=en 

in which I am collaborating with Martin to develop powerful new rules for integrating trigono- 
metric expressions. In particular for your next DNL, you might find of interest my last post in 
which I outline how recurrence equations are transformed into integration rules. 

Martin is certainly one of the world's most active and knowledgeable Derive users. Appar- 
ently he lives in Germany and is a member of DUG. Do you know his last name? 

Aloha, 

Albert 

Does anybody know, who “Martin” is. I found an email-address which could have been Mar- 
tin’s one: (clicliclic@freenet.de ). It would be great if we could find out who is the „secret“ 
DERIVE user? (The email address is not valid any more!) 


David Sjostrand from Sweden wrote: 

Dear Josef, 

after only having had a short glance at your and Carl's paper in DNL #80 I couldn't resist 
investigating some more sequences with Excel. (I shouldn't have done this. I have so many 
things I must do.) 

It seems that the sequence P(n)=m*P(n-3) + (l-m)*P(n-2) is convergent for all m and all ini- 
tial conditions. I attachjosef_carl.xls. 

David sent a proof for DNL#80, page 22 in his next mail and we started an exciting exchange 
of mails. The results of this correspondence will form an article in DNL#82, Josef 
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Notes by R.N. (Roger) Folsom 

Monterey, California, USA 
rnfolsom@redshift.com . 


#1: CaseMode := Sensitive 

As you will see below, I am trying to formally prove the circumstances under which a particu- 
lar inequality is true — and Derive's Simplify or Solve command actually says that it is true. 
(For now, all I need is sufficient conditions.) 

But I have not been able to figure out the Derive6.1 command or commands that will do that. 

My sources are Bernhard Kutzler and Vlasta Kolol-Voljc's "Introduction to Derive 6, Ad- 
vanced Mathematics for Your PC" (hereafter Derive6); and Albert and Joan Rich, Theresa 
Shelby, and David Stoutemeyer's "Derive User Manual Version 6" (hereafter Derive3). 

Derive6's equation solver (page 57), when used to determine whether a relationship is True 
or False, doesn't apply to what I am trying to do because it doesn't allow conditions to be 
imposed on the equation's (or, in my case, the inequality's) variables. 

Derive6's discussion of the IF tool (on pages 166 and 227) didn't help because I didn't un- 
derstand page 166 (I don't know any Analytic Geometry), and although I think that I did un- 
derstand page 227 I couldn't figure out how to apply it to my problem. 

Derive3's discussion of the IF tool (pages 309-310), made sense, but again I couldn't figure 
out how to apply it to my problem. 

My understanding is that the IF tool has four fields: a test condition, a "then" clause of re- 
sults if the test condition is satisfied, an "else" clause of results if the test condition is NOT 
satisfied, and an "unknown" clause. 

Derive3's opening IF example (page 309) is a payroll problem: IF(h <= 40, 1 0h, 400 + 15(h- 
40)) where h is hours worked and 1 0h and 15(h-40) are worker payment amounts. Unfortu- 
nately, the problem does not include the unknown clause, which if not included but is needed 
turns out to be the entire IF(...) tool. An example that does include an unknown clause would 
have been helpful. 

What I want is a tool, or set of tools, that let me state restrictions on my inequality's variables, 
then state the inequality itself, and returns either "true" or "false", or ? (or the same object) if 
the restrictions are too weak to determine whether the inequality is true or false. 


The rjn variables below are annual (or some other time period) rates of growth, such as in- 
terest rates, inflation and deflation rates, real output growth or decline rates, etc. Rates of 
growth rjn can be either positive, zero, or negative. 

The number N is the number of time periods being considered. Sometimes "total" rates of 
growth over N time periods are calculated by summing the N periodic rates of change rjn (n 

= 1,2 N), which is less accurate than multiplying the N corresponding (1+rjn) factors. 

Summing the rjn can give results that are either larger or smaller than the more accurate 
multiplication of the (1 + rjn). 
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However, a common situation is that over N time periods, more rates of growth rjn are posi- 
tive than negative. If that difference is sufficiently large, summing the N periodic rates of 
change rjn gives a lower total rate of growth over N periods than does multiplyng the corre- 
sponding (1 + rjn) numbers, so that summation understates the true total rate of change 
over N periods. 


Mathematically, that "common situation" inequality can be written in three basic ways. See 
statements #2, #3, and #4, below. 


SATISFYING THE "COMMON SITUATION" INEQUALITY; A SUFFICIENT CONDITION: All 
rjn > 0 

If ALL of the rjn are positive or zero (rjn > 0), the "common situation" inequality holds, be- 
cause the left side's terms are smaller than the right side's factors. For an example, let N = 
3, and see statements #5, #6, #7, and #8. Note that in #7 — a rearrangement of #2 with N = 
3 — the left side is -1 , and ALL of the variables on the right side are positive. 

Therefore, a SUFFICIENT condition for the "usual condition" inequality is simply that all all 
rjn >= 0. 


N N 

#2: I r < fl (1 + r ) 

n=l n n=l n 


#3: 


0 < 


N 

n 


(1 + r ) 
n > 


N 

I r 
n=l n 


#4: 


-1 < 


N 

n 

^n=l 


(1 + r ) 
n > 


1 


N 

I r 
n=l n 


For example, if N = 3: 


# 5 : 


-1 < 


3 

n 

<n=l 


Cl + r ) 
n j 


1 


3 

I r 
n=l n 


Statement #5 can be written as 

-1 < (r + l).(r + l).(r + 1) - 1 - (r + r + r ) 

#6 : 1 2 3 3 2 1 

In statement #6, multiplying out the (rjn + 1) factors would eliminate the -1 - (r|3 + rj2 + rjl) 
expression. Doing that (by using Derive's Simplify-Factor command) gives statement #7. 

-1 < r -(r -(r +l)+r)+r-r 

#7: 3 2 1 1 12 

or 

-1 < r ■ r ■ r + r ■ r + r ■ r + r ■ r 
#S: 123 23 13 12 


Note: I have done the same calculations above using N = 5. 
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Given the sufficient condition that all rjn > 0, statements #7 and #8 are obviously true be- 
cause the left side is negative, and every term on the right side is positive. 

And since statements #7 and #8 are merely a rearrangement of #5, which is an example of 
statements #4, #3, and #2, statement #7 demonstrates that if all rjn > 0, the "usual situation" 
inequality holds, so that summing N periodic rates of change rjn (instead of multiplying cor- 
responding 1 + rjn numbers) understates the true total rate of change over N periods. 

HOWEVER, I need to learn how to persuade Derive6.1 to do a formal proof (results either 
true, false, or unknown) of the preceding paragraph's results, because my next step is to 
experiment with less restrictive conditions, such as all rjn >= -1 . 


TO DO: Figure out a formal proof — for any statement #5 through #8 — that when #10 (be- 
low) holds, comes up with "true" rather than the wrong answer of "false." 

For example, Simplify-Basic of either IF(#10, #7) or IF(#7, #10) returns the same (written out) 
IF statement, presumably meaning that Simplify could not determine whether the inequality 
was true or false, despite the "all rjn >= 0" variable conditions. 

But Solve-Expression of either IF(#10, #7) or IF(#7, #10) concludes that the constrained ine- 
quality is false, when I know that it is true (from the analysis that follows #8). 

I clearly am not using IF properly. It may not be intended for what I am trying to do. In that 
case, I need to know what else to do. 


These are two Strange Attractors. I am quite sure that I was the 
first person who saw them on a PC screen. Josef 
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I received a mail from Robert Setif, a DUG Member from the very first times of its existence: 

Dear Josef, 

(in bad English) 

I am always alive and I spend my time to program often deranged counting with different mathemati- 
cal softwares (Derive, Mathematica, Maple, Mupad, Xcas, python, pari-gp). Regarding these soft- 
wares, do you know others, free preferably? 

I send you 3 files which calculate decimals of PI. You can make what you want.of them.I I am unable 
to ameliorate them on billing(displaying, printing ). In fact, the billing of long integers or floats with a 
lot of decimales by edges(slices) of 3 or 5 or 10 digits.lt would be more practical for checking. Could 
you fix me on this subject? 

I should like if there is a function that for instance M time_elapsed(pi_machin( 10001)) for calculating 
the time for that calcul. Thank you very much. I congratulate you for your work. 

Best regards. 

(in good French) 

Je suis toujours en vie et je passe mon temps a programmer des calculs souvent delirants avec diffe- 
rents logiciels mathematiques (Derive, Mathematica, Maple, Mupad, Maxima, Xcas, python, pari-gp). 
A propos de ces logiciels, en connaissez-vous d'autres, gratuits de preference? 

Je vous envoie 3 fichiers qui calculent des decimales de PI. Vous pouvez en faire ce que vous voulez. 
Je suis incapable de les ameliorer sur faffichage. En effet, faffichage de longs entiers ou de flottants 
avec beaucoup de decimales par tranches de 3 ou 5 ou 10 chiffres serait plus pratique pour la verifica- 
tion. Pourriez-vous me depanner a ce sujet? 

J’aimerais aussi une fonction comme par exemple M time_elapsed(pi_machin( 10001)) qui donnerait le 
temps necessaire pour ce calcul. Merci infiniment. 

Je vous felicite pour votre travail. Bravo ! 


These are Robert’s files which were attached to his mail : 


#1: InputMode := Word 


appin(x, nd) := 

Prog 

Precision := Approximate 
PrecisionDigits := nd + 5 
NotationDigits := nd 

nit := CEILING(nd/ (2 • L0G(l/x , 10))) + 2 
DISPLAY(nit) 

#2 : t := l/(2*nit - 1) 

k := 2«nit - 3 
Loop 

t := 1/k - t-x-x 
k := k - 2 
If k < 1 exit 
RETURN t • x 


#3: 


r 

appin 




1 

239 


x 

60 

J 


15 


#4: 


0 . 00418407600207472386453821495928545274104806530763195082701961 
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appi(x, nd) := 

Prog 

Precision := Approximate 
PrecisionDigits := nd + 5 
NotationDigits := nd 

nit := CEILING(nd/ (2 • LOG(l/x , 10))) + 2 
#B: t := 1/(2- nit - 1) 

k := 2- nit - 3 
Loop 

t := 1/k - t-x-x 
k := k - 2 
If k < 1 exit 
RETURN t • x 


# 6 : 


appi 


L 239 



#7 : 0 . 004184076002074723864B382149B928B4B274104806B3076319B082701961 


# 8 : 

#9: 



' 1 > 

r 1 i 


' i 

pig(nd) := 12«appi 

, nd 

+ 4*appi , nd 

+ 4*appi 

, nd 


l 4 J 

l 20 J 


l 198 B ) 


pig(101) 


#10 : 3 . 141B926B3B89793238462643383279B0288419716939937B10B820974944B9230781640628620899862- 

80 3482 B 3421170679 


#11: pig (1001) 

#12 : 3 . 141B926B3B89793238462643383279B0288419716939937B10B820974944B9230781640628620899862- 

803482B342117067982148086513282306647093844609BB0B822317253B94081284811174B02841027~ 
01938B2110BBB96446229489B493038196442881097B66B93344612847B64823378678316B271201909~ 
14B648B669234603486104B4326648213393607260249141273724B8700660631BB88174881B2092096~ 
28292 B409171B 364367892 B90 360011330 B 30 B48820466B213841469B1941B116094330B727036B7B9B- 
919B 309218611738193261179310 B118B48074462 379962749B67 3B188B7B2724891227938183011949- 
1298336733624406B664308602139494639B224737190702179860943702770B392171762931767B238- 
467481846766940B132000B6812714B263B6082778B7713427577896091736371787214684409012249- 
B3430146B49B8B3710B07922796892B8923B420199B61121290219608640344181B9813629774771309- 
960B187072113499999983729780499S10B9731732816096318B9B0244B94BB3469083026425223082B- 
334468B0352619311881710100031378387B2886B87533208381420617177669147303B982B34904287- 
BB4687311B9B628638823B3787B937B19B77818B7780B3217122680661300192787661119B909216420- 
1989 


#13: pig(2001) 


The second file was pig.dfw which was pretty the same as pigl .mth above including a com- 
ment that Robert had calculated n up to 5000 decimal places and needed approximatey 12 
minutes. 

I wrote back and asked for some information about this algorithm - and I had some ideas 
how to answer his request applying string manipulations and the time counter which was 
presented by Rene Hugelshofer some time ago. Using the DNL-index I found out that this 
was in DNL#55. 

But read first Robert’s information about this algorithm for calculating n with arbitrary accu- 
racy: 
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Dear Josef, 

Explications about functions appin et appi of my approx of PI. 

Of Aide for free software XCAS , aid written by Renee de GRAEVE Universite J Fourier 
GRENOBLE FRANCE. This AIDE is enormous. I have not the courage for translate it for you. You 
may find it on the site 

www-f ourier . u j f-grenoble . f r/~parisse/ giac_f r . html 

Bernard Parisse author of XCAS Universite J Fourier GRENOBLE FRANCE 
Please excuse that shambles of French-English. 

Atan(x)=x-x A 3/3+x A 5/5. . . PI/4=4 atan(l/5) - atan(l/239) 

Pour calculer la somme de n termes de la serie on va utiliser la methode de Horner pour faire le moins 
possibles de multiplications, on a : 

arctan(x)=x(l-x 2 (l/3-x 2 (l/5-x 2 (....-x 2 (l/2/i-l))))) 

L’utilisateur doit rentrer la valeur a de x et le nombre n de termes de la serie. 

gregory (a, n) :={ 
local t , k; 
t : = 1 / ( 2 *n-l ) ; 

for (k:=2*n-3; k>0; k:=k-2) { 

t :=l/k-a A 2 *t; 

} 

return (a*t) ; 

} ; 

On tape : 

16 * gregory (1/5,18) -4 * gregory (1/23 9, 6) 

On obtient : 

3.14159265359 

On tape : 

evalf (16*gregory (1/5,42) -4 * gregory (1/239,20) ) 

On obtient avec DIGITS :=60 

3.14159265358979323846264338327950288419716939937510582097494 

Soit R n (x ) le reste de la serie X&=i°° {-\) k+l x 2k [ !2k-\ : R n {x)=^ n+ ™(-\) k+x x 2k - x l2k- 1. 

On sait que |i?„(x)|<|x| 2 " +1 /2«+l Pour avoir |i?„(l/5)|=|i?„(0.2)|<10 -61 il faut que :2 2 " +1 <(2«+l)10 2 "“ 60 et 
comme 2'0^ 10 3 cela donne si on suppose 2«+l>10 : 

10 (6« + 3)/io <10 2«-59 soit 593<1 4„ soit 42 on verifie pour «=42 on a (l/5) 83 /85<2.56e-62 
On peut aussi ecrire si on suppose que 2«+ 1 > 1 0 : |x| 2,,+ 1 /2«+ 1 <x 2 " + 1 / 1 0< 1 0" 6 1 done on va choisir 
x 2 ” +1 <1CT 60 soit (2n+l)/ogl0(x)<-60 ou encore w>((-60)//ogl0(x)-l)/2 Pourx=l/5 on a 
«>42. 4202967422 et comme 2«+l>40 on peut ameliorer la majoration |x| 2 " + 1 /2h+ 1 <x 2 ” + 1 /40< 1 0 6 1 ce 
qui donne «>((-60+/ogl 0(4))//ogl 0( l/5)-l )/2=41 .989620 1 841 done on prend //=42 pour x=l/239 on 
a n>l 1.66 done on prend «=12 et on verifie : (l/239) 25 /25=l .3871 1499837^—6 1 . 

On choisit 62 Chif f res pour DIGITS et on tape : 

evalf (16*gregory(l/5, 42) -4*gregory (1/239, 12) ) 

On obtient : 

3.1415926535897932384626433832795028841971693993751058209749446 
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On tape : 

evalf (pi 

On obtient : 

3.1415926535897932384626433832795028841971693993751058209749446 

Remarque 

Avec cette methode John Machin calcula 100 decimates de n en 1706. 

Then I asked Robert how the result should look like and presented some alternatives. His 
answer was as follows: 

Dear Josef, 

It would be like 3.14159 26535 89793 23846 26433 
Concerning the calculation time, 3 possibilities: 

1. calculate pig(1001) in the background - without showing the result, 

2. see only the calculation time as result 

3. see the calculation time stored in a global variable 

These 3 solutions are very nice, but I can wait patiently... (I was born in 1929 during the crysis!) 
Best regards. 

Robert 


First of all my congratulations to Robert. I really hope that I when I will reach 82 years of life 
time will also be interested in presenting some 1000 digits of tt!! 


This is what I found out finally (using Rene Hugelshofer’s time couting trick from DNL#55): 

#1: InputMode := Word 

Global Variables: 

#2: [calc_time :=, pi_res :=] 

appi(x, nd) := 

Prog 

Precision := Approximate 
Preci sionDigits := nd + 5 
NotationDigits := nd 

nit := CEILING(nd/(2.L0G(l/x, 10))) + 2 
#3: t := 1/(2- nit - 1) 

k := 2. nit - 3 
Loop 

t := 1/k - t • x • x 
k := k - 2 
If k < 1 exit 
RETURN t • x 


#4: pig(nd) := 32«appi 


* , nd 


l 10 


r i 


- 4*appi 


l 239 


■ , nd 


r 1 


16«appi 


L 515 



My first attempt which did not satisfy Robert s expectations: 


#5: pig_time(100 , 20) 


# 6 : 


T 3.141592653589793238 
46264338327950288419 
71693993751058209749 
44592307816406286208 
99862803482534211706 
_L 7 


0.05 
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Output in columns 

pig_t_col (nd p digs, t, t_, time, res, out) := 

Prog 

t_ := RANDOM (0) 

res := STRING(32 -appi (1/10 , nd) - 4-appi (1/239 , nd) - 16- appi (1/515 , nd)) 
time := APPROX((RANDOM(0) - 0/100) 
calc_time := time 
If DIM(res) < digs 
RETURN [res, time] 

#7 : out := [ ] 

Loop 

out := APPEND(out, [res|[l, .... digs]]) 

res ;= res|[digs + 1 DIM(res)] 

If DIM(res) < digs 
If t = 1 

RETURN [ [APPEND(out , [res])] 1 , APPEND( M Calc time = ", STRING(time))] 
RETURN [ [APPEND(out , [res])] 1 ] 


Output in columns to "digs" digits without calculation time, 
Calculation time is stored in global variable "calc_time" 
which can be given after the calculation (#10) 


#8: 


#9: 


pig_t_col (501, 40) 

T 3.14159265358979323846264338327950288419 
7169399375105820974944592307816406286208 
9986280348253421170679821480865132823066 
4709384460955058223172535940812848111745 
0284102701938521105559644622948954930381 
9644288109756659334461284756482337867831 
6527120190914564856692346034861045432664 
8213393607260249141273724587006606315588 
1748815209209628292540917153643678925903 
6001133053054882046652138414695194151160 
9433057270365759591953092186117381932611 
7931051185480744623799627495673518857527 
2489122793818301194912 


#10: calc_time = 1.43 


Output in colums with 1 as third parameter, calculation time is given in the output 


# 11 : 


pig_t_col (201, 40, 1) 

T 3.14159265358979323846264338327950288419 ' 
7169399375105820974944592307816406286208 
9986280348253421170679821480865132823066 
4709384460955058223172535940812848111745 
0284102701938521105559644622948954930381 
1 96 


Calc time = 0.17 


# 12 : 
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Same procedure but output in rows (and "digs" slices) with default digs = 5: 

pig_t_row(nd , digs := 5 , t, t_, time, res, out) := 

Prog 

t_ := RANDOM(O) 

res := STRING(32 *appi (1/10 , nd) - 4-appi (1/239 , nd) - 16 -appi (1/515 , nd)) 
time := APPROX((RANDOM(0) - 0/100) 
calc_time := time 
If DIM(res) < digs 
RETURN [res, time] 

#13: out := "3." 

res := REST(REST(res)) 

Loop 

out := APPEND(out, res|[l digs], " ") 

res := resj[digs + 1 DIM(res)] 

If DIM(res) < digs 
If t = 1 

RETURN APPEND(out , res, ", Calc time = ", STRING(time)) 

RETURN APPEND(out , res) 

#14: pig_t_row(101, 5) 

#15: 3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37510 58209 74944 59230 78164 


06286 20899 86280 34825 34211 70679 


#16: calc_time =0.04 

#17: pig_t_row(101, 5, 1) 

#18: 3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37510 58209 74944 59230 78164 
06286 20899 86280 34825 34211 70679, Calc time = 0.05 


#19: pig_t_row(501, 20) 

#20: 3.14159265358979323846 26433832795028841971 69399375105820974944 59230781640628620899 

86280348253421170679 8214808651328230664 7 09384460955058223172 53594081284811174502 

84102701938521105559 64462294895493038196 44288109756659334461 28475648233786783165 

27120190914564856692 34603486104543266482 13393607260249141273 72458700660631558817 

48815209209628292540 91715364367892590360 01133053054882046652 13841469519415116094 

33057270365759591953 09218611738193261179 31051185480744623799 62749567351885752724 

89122793818301194912 
#21: calc_time =1.46 

pig_time_only(nd , digs, t_, time, res, out) := 

Prog 

t_ := RANDOM(O) 

#22: pi_res := 32 -appi (1/10 , nd) - 4-appi (1/239 , nd) - 16 *appi (1/515 , nd) 

time := APPROX((RANDOM(0) - 0/100) 

RETURN time 

The last function gives back the calculation time only, 
the resulting number is stored as "pi_res" 

#23: pig_time_only(501) 

#24: flEW 

#25: pi_res = 

3.141592653589793238462643383279502884197169399375105820974944592307816406286208998- 

62803482534211706798214808651328230664709384460955058223172535940812848111745028410- 

27019385211055596446229489549303819644288109756659334461284756482337867831652712019- 

09145648566923460348610454326648213393607260249141273724587006606315588174881520920- 

96282925409171536436789259036001133053054882046652138414695194151160943305727036575- 

95919530921861173819326117931051185480744623799627495673518857527248912279381830119- 

4912 


#26: pig_time_only(2001) 

#27 : 60.36 


Robert’s answer: Thank you very much for your work. It is great (magnifique en frangais). 
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Mathematical Physics Electronic Journal , University of Barcelona 
Probability Surveys 

Publications de llnstitut Mathematique, Mathematical Institute of the Serbian Academy of Arts and 
Sciences 

Rendiconti, University of Torino 

Revista Colombiana de Estadlstica (Colombian Journal of Statistics) 

Revista Colombiana de Matematicas, University of Bogota 
Seminaires & Congres, proceedings of mathematical meetings 

Seminaire Lotharingien de Combinatoire, Universities of Bayreuth, Erlangen and Strasbourg 

Siberian Electronic Mathematical Reports, Sobolev Institute of Mathematics 

Surveys in Mathematics and its Applications, University Constantin Brancusi of Targu-Jiu 

Symmetry, Integrability and Geometry: Methods and Applications (SIGMA), Ukraine 

The Journal of Nonlinear Science and its Applications, Iran 

Theory and Applications of Categories, Sackville, Canada 

TICMI, Tbilisi International Centre of Mathematics and Informatics 

Universitatis Iagellonicae Acta Mathematica, University of Krakov, Poland 

Journal of Analysis and its Applications, Lemgo, Germany 

Mathematical Proceedings of the Royal Irish Academy 

Memoirs on Differential Equations and Mathematical Physics, Tbilisi 

The New York Journal of Mathematics 

Novi Sad Journal of Mathematics 


The EMIS-website http://www.emis.de/elibrn/ is hosting also a rich collection of 

• Conference Proceedings 

• Monographs and Lecture Notes 

• Classical Works, Selecta and Opera Omnia (among others Riemann and Hamilton) 

• Software and Other Special Electronic Resources (among others: Mathematical English 
Usage - A Dictionary) 


